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A few topics concerning the early chemical evolution of the Milky Way are critically discussed. 
In particular, it is argued that: 1) Observed abundance patterns of extremely metal poor stars (of 
Pop. II) do not constrain the mass range of the first generation (Pop. Ill) stars; the latter may 
well be normal massive stars (10-50 M ) or very massive ones (140-1000 M Q ) or a combination 
of the two classes. 2) The discrepancy between primordial Li abundance (after WMAP) and the 
observed "Spite plateau" cannot be due to astration by a generation of massive Pop. Ill stars, as 
recently suggested, unless if such stars eject negligible amounts of metals. 3) The observed halo 
metallicity disribution may well be understood in the framework of hierarchical galaxy formation, 
as shown here with a simple semi-analytical model. 4) Formation of the Milky Way's halo from 
a myriad of smaller sub-haloes may have important implications for our understanding of the 
abundance patterns of r-elements, the origin of which remains still unclear. 
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1. Introduction 

Galactic chemical evolution studies provide a useful framework for the interpretation of the 
everincreasing wealth of observational data concerning abundances in stars of the Milky Way. In 
particular, studies of abundance ratios in the oldest stars of the Milky Way (Pop. II stars, with 
metallicities down to [Fe/H]=log(Fe/H)-log(Fe/H) =-4.) allow one to put interesting constraints 
on the nature of the first stars that enriched the interstellar medium with metals. Ultimately, one 
may hope to find a distinctive imprint of the very first generation of massive stars (Pop. Ill), those 
that were born out of primordial material containing H and He only (as well as trace amounts of 
7 Li, resulting from the Big Bang). This pattern should then be visible on the surfaces of the next 
generation stars that were born from the enriched gas and are still around today, i.e. the oldest and 
most metal poor Pop. II stars. 

2. Halo star composition: any signatures of massive Pop. Ill stars ? 

Star formation theory and numerical simulations suggest that the first stars were very massive 
(e.g. Nakamura and Umemura 2002, Glover 2004 and references therein). In view of the many 
important implications of those ideas (e.g. for re-ionization of the Universe) it is interesting to 
investigate whether the observed abundance patterns of old stars reveal any trace of a Pop III 
stellar generation with peculiar properties. In fact, it has been suggested that a stellar initial mass 
function (IMF) rich in massive stars (top-heavy IMF) can help explain (1) the abundance ratios 
in the Fe-peak region, and (2) the discrepancy between the expected primordial Li value (after 
the determination of the cosmic baryonic density by WMAP) and the one observed in the "Spite 
plateau". We argue below that the answers are: (1) may be, and (2) no. 

2.1 Abundance patterns in EMP stars: no clear signature of the culprits ! 

The recent data on Extremely metal Poor (EMP) stars confirm some previously found tenden- 
cies in the abundance patterns of heavy elements, and in particular: 

- The constancy, down to the lowest observed metallicities, of the a/Fe ratio (where a stands 
for alpha elements like O, Mg, Si, Ca); the high value of that ratio (<~3 times solar) is attributed to 
the late contribution of Fe by SNIa, with long lived progenitors. 

- The decreasing trend of Mn/Fe and Cu/Fe with decreasing metallicity, in quantitative agree- 
ment with theoretical expectations (see Fig. 1). 

However, more interesting are the cases where observations are at odds with theoretical ex- 
pectations. These are the cases of Cr, Co and Zn, as seen in Fig. 1, where the results of a full-scale 
chemical evolution model of the Milky Way (originally presented in Goswami and Prantzos 2000) 
are displayed. The adopted stellar yields are from non-rotating Core Collapse supernovae (CCSN) 
in the "normal" mass range of 12 to ~40 M , exploding with "canonical" energies of E =10 51 erg 
and with progenitor metallicities from Z=0 to Z=Z . 

The reasons for those discrepancies probably lie in our poor understanding of the explosion 
mechanism and of the nature of the early CCSN. A large amount of work was devoted in recent 
years to understanding the peculiar abundance patterns among Fe-peak elements observed in EMP 
stars (e.g. Heger and Woosley 2002, Maeda and Nomoto 2003, Chieffi and Limongi 2004; for 
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Figure 1: Abundance ratios [X/Fe] as a function of metallicity [Fe/H] in stars of the Milky Way; small data 
points are from various sources, while the large data points at low metallicity are from the VLT survey of 
Cayrel et al. (2004). Results of standard chemical evolution models, performed with two sets of metallicity 
dependent massive star yields (Woosley and Weaver 1995, solid curves ; Chieffi and Limongi 2004, dashed 
curves) are also displayed. Yields for SNIa are adopted from Iwamoto et al. (1999) and for intermediate 
mass stars from van den Hoek and Gronewegen (1997). While the behaviour of alpha elements and Mn is 
correctly reproduced (at least qualitatively), there are large discrepancies in the cases of Cr, Co and Zn. 



reviews see Nomoto et al. 2005, Cayrel 2006 and references therein). In those works it was ex- 
plored whether some specific property of the Pop. Ill stars (characteristic mass, rotation, explosion 
energy etc.) was quite different from the corresponding one of their more metal rich counterparts 
and produced a different nucleosynthetic pattern in the ejecta. Two main classes of solutions (not 
mutually exclusive) appear to arise from those works: 

1) Normal massive stars (10-50 Mq), exploding as core collapse SN (CCSN) and leaving 
behind a neutron star or a black hole. It was shown (essentially by the Tokyo group, see Nomoto 
et al. 2005 for a recent review) that higher energies than the canonical one of E = 10 51 ergs 
combined with asphericity of the explosion can help to improve the situation concerning the Zn/Fe 
and Co/Fe ratios observed in EMPs. Material ejected along the rotation axis (in the form of jets) 
has high entropy and is found to be enriched in products of a-rich freeze-out (Zn and Co), as 
well as Sc (which is generically undeproduced in spherical models); this kind of models seems 
at present the most promising, but this is not quite unexpected (since they have at least one more 
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Figure 2: Nucleosynthesis yields of a 25 M© star, compared to observations of extremely metal-poor stars 
(vertical error bars, from Cayrel et al. 2004). Left: The star explodes with a canonical energy of 1.5 10 51 
erg; Right: The star explodes as a hypernova, with an energy of 10 10 51 erg (both fi gures from Nomoto et al. 
2005). Abundance ratios in the Fe-peak are much better reproduced in the second case. 



degree of freedom w.r.t. spherical models). Such properties (high energies, asphericity) are indeed 
observed in the local Universe for hypernovae, a class of energetic CCSN. However, hypernovae 
are not so frequent today as to affect the late chemical evolution of galaxies (normal CCSN plus 
SNIa can easily account for e.g. the solar abundance pattern). If hypernovae indeed affected 
the early chemical evolution, they muct have been much more abundant in Pop. Ill than today. 
No satisfactory explanation for such a large hypernova fraction in the early Universe exists up to 
now. It seems hard to conceive that the typical energy of the explosion was substantially higher in 
Pop. Ill stars (all other parameters, except metallicity, being kept the same). On the other hand, 
rotation effects may indeed be stronger in low metallicity environments, because of smaller angular 
momentum losses due to lower mass loss rates (e.g. Meynet et al. 2006). In any case, none of the 
models explored sofar appears to account for the early trend of Cr/Fe, which remains a mystery at 
present. 

2) Very massive stars (VMS), above 100 M Q . Such stars are thought to: collapse to black holes, 
if M<140 M ; explode during oxygen burning as pair-implosion SN (PISN), if 140<M/M Q <300; 
and again collapse without ejecting metals if M>300 M (see Heger and Woosley 2002 and refer- 
ences therein). In the case of PISN, it was shown that they do not produce enough Zn to account for 
the oberved high Zn/Fe ratio in EMP stars (Heger and Woosley 2002, Umeda and Nomoto 2005) 
and they were thus excluded as major contributors to the early chemical evolution of the MW. This 
shortcoming cast doubt as to the existence of a Pop. Ill generation composed exclusively of VMS. 
However, in a recent work Ohkubo et al. (2006) explore the conditions under which rotating stars 
of M=500 and 1000 M do explode and eject metals (although a massive black hole is ultimately 
produced in both cases). They find explosions for a certain region of their parameter space and 
substantial production of Zn as a generic feature (see Fig. 3b). Even if each class of PISN or 
500-1000 M Q stars cannot reproduce the observed EMP abundance pattern alone, it is obvious that 
considering the full 140-1000 M range (folded with an appropriate IMF) will account for such a 
pattern. Thus, VMS cannot be excluded at present as candidates for Pop. Ill stars, at least not on 
"chemical" grounds. 
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Figure 3: Nucleosynthesis yields of Very Massive Stars (VMS, with mass >100 M©). Left: Stars in the 
140 to 260 M© range explode as Pair Instability Supernovae (PISN) and do not produce enough Zn (from 
Heger and Woosley (2002). Right: Stars in the 300-1000 M© range explode after core collapse and produce 
significant amounts of Zn (from Ohkubo et al. 2006). Abundance ratios in the Fe-peak are much better 
reproduced in the second case. A Pop. Ill population composed exclusively of VMS (in the 140-1000 M© 
range) and with an appropriate IMF is obviously compatible with observations of EMP stars. 

2.2 WMAP Li was (most probably) NOT astrated in massive Pop. Ill stars 

The recent WMAP results from the 3-year data analysis of the cosmic microwave background, 
combined with other cosmological measurements, allow to determine with great precision the pa- 
rameters characterising our observable Universe (see Lahav and Liddle 2006). In particular, the 
baryonic density is 1^=4.5 10~ 2 . For that density, calculations of Standard Big Bang Nucleosyn- 
thesis (SBBN) predict an amount of primordial deuterium that is fully compatible with observations 
of its abundance in remote gas clouds (at such early times that significant astration had no time to 
occur); this undoubtely constitutes a triumph for SBBN (e.g. A. Coc, this Conference). Unfortu- 
nately, the situation with Li ( 7 Li) is far less satisfactory. 

Indeed, the Li abundance of the "Spite plateau" (Li/H~ 1.5 10~ 1() =const. for halo stars, down 
to the lowest metallicities) is a factor of 2-3 lower than the WMAP+SBBN value. Barring sys- 
tematic errors, the conclusion is that primordial Li has been depleted, either before getting into the 
stars it is observed today, or during the lifetime of those stars. Two "depletion agents" have been 
proposed in the former case: decaying supersymmetric particles (Jedamzik 2006) and astration in a 
first generation of exclusively massive (mass range m*=10-40 M Q ) Pop. Ill stars (Piau et al. 2006). 
The latter idea, however, suffers from a serious flaw, since in that case the metallicity of the ISM 
(out of which the next stellar generation would form with depleted Li) would necessarily rise to 
levels much higher than those observed in EMP stars. This can be seen as follows: 

Assuming that the current halo stellar mass (M#=2 10 9 M Q ) was initially in the form of gas, a 
fraction / of which was astrated through massive stars, the resulting Li mass fraction is Xu=(l — 
f)Xu,p where Xuj> is the primordial Li abundance (assuming a return fraction R ~ 1 for the astrating 
and metal producing stars). Similarly, the resulting oxygen abundance would be Xo=niolMH, where 
the mass of oxygen mo = Nsn Yo is produced by a number of supernovae Nsn= /Mm / m* , each one 
with a typical oxygen yield of Yo (in M , to be discussed below). Putting all this together, one has: 

Xu_ = 1 O2g (X o /O.OO7)K/4OM ) 
Xn,p Yo 
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Figure 4: Left: Observations of Li in halo and disk stars of the Milky Way. The primordial Li value, obtained 
from the baryonic density of WMAP and calculations of SBBN, is indicated by a horizontal dashed line. 
The observed Li plateau at low metallicities depends sensitively on assumed stellar temperature, and differs 
from the WMAP value by factors 2-3. Right: Illustration of the problem encountered by the idea that the 
discrepancy is due to astration of Li by a Pop. Ill composed of exclusively massive stars, in the 10-100 
M Q range (Piau et al. 2006): such stars necessarily eject metals, either through their winds (e.g. Nitrogen, 
in case of rotating stars) or through the final supernova explosion (e.g. Oxygen). Mixing the ejecta with 
various proportions of primordial material would result in Li depletion (by a factor of ^2 in case of a 50-50 
mixture), but the resulting metallicity would be much higher than the one of EMP stars unless abnormally 
low oxygen yields were assumed (i.e. the curve in Fig. 4b parametrized with a yield io=0.001 M Q , see text). 



That relation appears in Fig. 4 (right panel) as a function of log(Xo/Xo ), with solar Xo Q =0.007. 
The four curves correspond to different assumptions about the typical oxygen yield of a massive 
star of Z=0, ranging from 0.001 to 1 M . Among the various calculations of massive star nu- 
cleosynthesis made sofar, only one set of models in the 30-40 M range (set A of Woosley and 
Weaver 1995, WW95) produces extremely low oxygen (and heavy element) yields for stars of 35 
and 40 M . This is due to the assumption of constant kinetic energy at infinity (=1.2 10 51 erg) 
adopted for all models of set A in WW95; other assumptions (sets B and C) lead to much larger 
oxygen yields and to abundance ratios in rough agreement with observations in EMP stars, at least 
for alpha-elements (e.g. Goswami and Prantzos 2000), while set A provides a rather poor match. 
Other calculations at Z=0, either for spherically symmetric stars (Chieffi and Limongi 2004, as- 
suming a fixed mass of 56 Ni ejected) or for bipolar explosions (Maeda and Nomoto 2003) lead to 
Yo > 1 M always; note that the latter models provide a better (albeit not perfect) match to observed 
Fe-peak element abundance ratios than "standard" (spherically symmetric) models, as dicussed in 
Sec. 2.1. 

Another way to eject astrated material by Z~0 stars is through stellar winds, which require 
rapidly rotating stars, since radiative pressure is innencient at such low metallicities; but rotating 
massive stars produce large amounts of nitrogen (which may in fact help explaining the observed 
primary-like nitrogen in EMP stars, e.g. Meynet et al. 2006), thus the problem of metal overpro- 
duction is not avoided in that case either. 

Astration in massive Pop. Ill stars cannot solve the Li discrepancy between the Spite plateau 
and WMAP+SBBN 1 : even a small Li depletion should be accompanied by excessive metal en- 

1 At least, not the 10-40 M, ;i stars suggested in Piau et al. 2006 (provided that current nucleosynthesis models for such 
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Figure 5: Left: Metallicity distribution of field halo stars from Ryan and Norris (1991, dots), and the 
ongoing research of Beers and collaborators (T. Beers, private communication, asterisks). The curve is a 
simple model with outflow rate equal to 8 times the star formation rate. Right: Stellar metallicity vs stellar 
mass for neary galaxies. Data and model are from Dekel and Woo (2003). The MW halo, with average 
metallicity [Fe/H]~-1.6 (see left panel) and estimated mass 2 10 9 M falls below that relationship. 

hancement. It appears now that depletion during the stellar evolution, within the stellar envelope, 
is the answer (see Korn et al. 2006 and references therein). 

3. The early MW and hierarchical galaxy formation 

In the previous sections, the early chemical evolution of the MW was discussed indepen- 
dently of the cosmological framework in which it took place. According to the currently dominant 
paradigm of hierarchical structure formation, the early phases of a galaxy's evolution are the most 
complex ones, as they involve multiple mergers of smaller sub-units. In the case of the Milky Way, 
interesting "chemical signatures" of that period should still be left around us today, in the form 
of the metallicity distribution (MD) of long-lived stars and in dispersion in abundance ratios. We 
discuss those two topics in the following. 

3.1 The halo metallicity distribution: outflow vs. subhalo merging 

The MD of Galactic halo field stars (HMD) is rather well known in the metallicity range - 
2.2<[Fe/H]<-0.8, while at lower metallicities its precise form has still to be established by ongoing 
surveys (Fig. 5a). Its overall shape is well fitted by a simple model of GCE as dn/dlogZ oc 
Z/y e~ z l y , where y is the yield of a stellar generation; this function has a maximum for Z = y. The 
HMD peaks at a metallicity [Fe/H]=-1.6 (or [0/H]=-l.l, assuming [O/Fe]~0.5 for halo stars), 
pointing to a low yield y=\l\3 of the corresponding value for the solar neighborhood. Such a low 
halo yield is "classicaly" (i.e. in the monolithic collapse scenario) interpreted as due to outflow 

stars are correct); 100 M stars collapsing to black holes would be better candidates (provided they eject a substantial 
fraction of their astrated mass). 
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[Fe/H] Mass of subhalo (M Q ) 



Figure 6: Left top and middle panels: Metallicity distribution (in lin and log scales, respectively) of the MW 
halo, assumed to be composed of a population of smaller units (sub-haloes). The individual MDs of the 
sub-haloes, from 10 6 M0to 2 10 8 MQare indicated in both panels (but clearly seen only in the middle), as 
well as the sum over all haloes (upper curves in both panels, compared to observations). Small sub-haloes 
contribute the largest fraction of the lowest metallicity stars (bottom left). Properties of the sub-haloes as a 
function of their mass. 

during halo formation (Hartwick 1976), at the large rate of 8 times the SFR (Prantzos 2003). How 
can it be understood in the framework of hierarchical merging ? 

It should be noted that the typical halo metallicity ([Fe/H]=-1.6) is substantially lower, by more 
than 0.5 dex, than the corresponding metallicities of nearby galaxies of similar mass {Mu a io ~2 
10 9 M ), as can be seen in Fig. 5b. That figure also displays the well known galaxian relationship 
between stellar mass and stellar metallicity; most probably, that relationship results from mass loss, 
which is more important in lower mass galaxies, since the hot supernova ejecta escape more easily 
their swallow potential well (e.g. Dekel and Woo 2003). 

Assuming that the MW halo has been assembled from sub-units similar to the low mass galax- 
ies of Fig. 5b, one may interpret the HMD as the sum of the MD of such low mass galaxies; it is 
assumed that each one of those galaxies evolved with an appropriate outflow rate and correspond- 
ing effective yield y(M) = (1 — R)/(l +k — R), where k(M) is the outflow rate in units of the SFR 
and R is the return mass fraction, depending on the adopted stellar IMF (see Prantzos 2003). In 
that case, one has: HMD(Z) = \IM Ha i JZ/y(M) e- z M M ^®(M)MdM, where 5>(M) is the mass 
function of the sub-units and y(M) is the effective yield of each sub-unit (obtained directly from 
Fig. 5b as y(M) = Z(M), i.e. smaller galaxies suffered heavier mass loss). 

The results of such a simple toy-model for the HMD appear in Fig. 6a. The HMD is extremely 
well reproduced, down to the lowest metallicities, assuming <I>(M) oc M~ 2 , such as those resulting 
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from recent high resolution numerical simulations for Milky Way sized dark haloes (Diemand et 
al. 2006, Salvadori et al. 2006). Low mass satellites (down to 10 6 M ) contribute most of the low 
metallicity stars of the MW halo, whereas the high metallicity stars originate in a couple of massive 
satellites with M~ 10 8 M (Fig. 6a, bottom). 

Some properties of the sub-haloes as a function of their mass appear also in Fig. 6 (right 
panels). The outflow rate, in units of the corresponding SFR, is k(M) = (1 — R){y~True/y(M) — 1), 
where yjme is the yield in the solar neighborhood (y t me = O.7Z , from the local G-dwarf MD). If 
the MW halo were formed in a a potential well as deep as those of comparable mass galaxies, then 
the large outflow rate required to justify the HMD (k=S) is puzzling; on the contrary, the HMD is 
readily understood if the MW halo is formed from a large number of smaller satellites, each one of 
them having suffered heavy mass loss according to the simple outflow model. 

A more physical, but much less "transparent", approach consists in deriving the full merger 
tree of the MW halo and (by using appropriate receipes for SFR and feedback for the sub-haloes) 
following the chemical evolution through merging/accretion with Monte-Carlo simulations. Recent 
studies (Tumlinson 2006, Salvadori et al. 2006) find good agreement with the observed HMD, but 
it is hard to find (in vue of the many model parameters) what is (are) the key factor(s) determining 
the final outcome. In any case, mass loss (occuring mostly through tidal stripping than through 
supernova feedback, e.g. Bekki and Chiba 2001) is crucial in shaping the HMD. 

3.2 Implications: abundance dispersion and origin of the r- elements 

The early chaotic history of the MW (according to the hierarchical merging scenario) bears 
naturally to the issue of dispersion in the abundance ratios of EMP stars. Simple-minded arguments 
suggest that dispersion should increase at low metallicities, under the explicit assumption that low 
metallicities correspond to such early times that complete mixing of SN ejecta with the interstellar 
medium is impossible (Argast et al. 2002, Karlsson and Gustafsson 2001). In order to actually 
see such a dispersion, variations in abundance ratios among SN of different masses should be 
sufficiently large (say, larger than typical observational statistical uncertainties of 0.1 dex). 

Several recent studies of EMP stars reveal very small scatter in the abundance ratios of all 
elements in the C to Fe-peak mass range (Carretta et al. 2002, Cayrel et al. 2004, Barklem et 
al. 2006). This could mean that i) mixing timescales are (much) shorter than typical chemical 
evolution timescales at metallicities down to [Fe/H]~-3.5 or ii) variations in abundance ratios of 
SN of different masses and energies (for the C to Fe-peak mass range) are sufficiently small, or a 
combination of (i) and (ii). Despite our current ignorance of the timescales of mixing and early 
chemical evolution, case (i), i.e. extremely efficient and rapid mixing of SN ejecta, cannot be the 
whole truth, since large variations in abundance ratios are observed in the case of r-elements. A 
solution to that may be that r-elements originate mostly in a restricted sub-class of CCSN (say 8-10 
M ), which produce a very large r/Fe ratio (e.g. Ishimaru et al. 2003, also these proceedings). 

The hierarchical merging scenario may complicate even more the already unclear situation of 
the r-element production. At present, the production site of r-nuclei is unknown, with (some class 
of) CCSN favoured relative to neutron star mergers (NSM). The reasons for the "disfavor" of NSM 
are explored and clearly presented in Argast et al. (2004): the combination of expected low event 
rates, long coalescence timescales and high individual yields leads to a relatively late injection of 
the NSM r-nuclei in the Galaxy (around [Fe/H]~-2.5, while Eu is observed already at [Fe/H]~-3 
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Figure 7: Illustration of how differently evolving sub-haloes may create a dispersion in the abundance ratio 
[X/Fe], provided element X is produced in two different nucleosynthesis sites, with different evolutionry 
timescales; this might explain observed dirpersion in r/Fe. Left: Halo 1 evolves rapidly, so that only the 
first (short-lived) site of element M operates, and it reaches a high metallicity (say [Fe/H]=-1.2); halo 2 (less 
massive) evolves on longer timescales, allowing for the second source of element M to increase the M/Fe 
ratio, but it reaches a smaller metallicity (say [Fe/H]=-2.6). Right: When those sub-haloes merge, dispersion 
as a fuction of metallicity is naturally obtained (starting at low [Fe/H]), even if the ISM in each one of them 
is completely homogeneized. 

and the r-component of Ba at even lower metallicities), and produces a large scatter of r/Fe ratio 
relatively late (at [Fe/H]~-2, when observations show that such large scatter occurs only at lower 
metallicities, around [Fe/H]~-3, Fig. 4 in Argast et al. 2004). 

The conclusions of the Argast et al. (2004) study depend a lot on the adopted assumptions 
about the NSM properties, as the authors recognize. They also depend on their model of the early 
Galaxy evolution, in which star formation proceeds stochastically in various places, but at the 
same average rate everywhere, i.e. there is in reality one single system evolving. The situation is 
obviously different in a (presumably more realistic) system composed of many different sub-units, 
each one evolving at its own rate and affected in a different way by mass loss. 

An illustration of such a situation is provided in Fig. 7: Halo 1 (HI, massive) evolves rapidly 
(100 Myr) and reaches a high metallicity ([Fe/H]~-1), while halo 2 (H2, low mass) evolves slowly 
(~1 Gyr) and reaches a lower final metallicity ([Fe/H]~-2.6). Because of short timescales, NSM do 
not enrich HI with r-elements, so its r/Fe ratio remains [r/Fe]=const=0. In H2, NSM start increasing 
the r/Fe ratio only after about 100 Myr (assumed as typical coalescence timescale of NSM), up 
to [r/Fe]~l (before 100 Myr, CCSN produce [r/Fe]=0). Note that efficient and rapid mixing of 
CCSN and NSM products in the gas of both haloes is always assumed (to ensure no dispersion 
in abundance ratios within each halo. When the two haloes merge together, two distinct histories 
appear in the r/Fe vs Fe/H plane (Fig. 7 right). By considering a large number of intermediate mass 
haloes, the plane can be naturally filled at [Fe/H]>-3. A large dispersion can then be naturally 
recovered and NSM can be important sources of r-nuclei (along with CCSN), while at the same 
time dispersion for other elements (exclusively produced in CCSN) remains small, in agreement 
with observations. 

In summary, because of the chaotic early history of the MW (formed by a myriad of smaller 
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sub-haloes with widely differing SF histories), it cannot be excluded at present that both CCSN 

and NSM have contributed to its enrichment with r-elements. 
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